The Charge Form Factor of the Neutron at Low Momentum Transfer from the 

2 H(e, e'n)p Reaction 
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We report new measurements of the neutron charge form factor at low momentum transfer using 
quasielastic electrodisintegration of the deuteron. Longitudinally polarized electrons at an energy 
of 850 MeV were scattered from an isotopically pure, highly polarized deuterium gas target. The 
scattered electrons and coincident neutrons were measured by the Bates Large Acceptance Spec- 
trometer Toroid (BLAST) detector. The neutron form factor ratio G%/Ctm was extracted from 
the beam-target vector asymmetry A^d at four-momentum transfers Q 2 = 0.14, 0.20, 0.29 and 0.42 
(GeV/c) 2 . 
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The neutron is composed of charged constituents, 
whose net distribution is described by the charge (or elec- 
tric) form factor G^. Differences in the up and down 
quark distributions produce a nonuniform distribution 
of the net charge The neutron electric form factor 
G% exhibits a maximum in the region of Q 2 « 0.1 — 0.5 
(GeV/c) 2 ; when Fourier-transformed this corresponds to 
a positively charged core and a concentration of nega- 
tive charge at intermediate to large distances of « 1 fm, 
commonly associated with a meson cloud surrounding 
the nucleon. Models which emphasize the role of the 
meson cloud have been successful in explaining impor- 
tant aspects of nucleon structure @, d, 0, S 0] • Precise 
knowledge of G^ is essential for the description of electro- 
magnetic structure of nuclei, and for the interpretation 
of parity violating electron scattering experiments to de- 
termine the strangeness content of the nucleon. Further, 
it is anticipated that exact ab initio QCD calculations 
of using lattice techniques will eventually be possi- 
ble 0j. 

In the absence of a free neutron target, determina- 
tions of G% at finite Q 2 are typically carried out using 
quasielastic electron scattering from deuterium or 3 He 
targets. Despite the small value of G% compared to the 



neutron magnetic form factor Gjjf, it can be obtained 
with high precision from double-polarization observables 
based on the interference of G^ with . With the avail- 
ability of high-duty-factor polarized electron beams over 
the last decade, experiments have employed recoil po- 
larimeters, and targets of polarized 2 H and 3 He to per- 
form precision measurements of G^ using polarization 
techniques with inherently small systematic uncertain- 
ties 0. The slope of G n E {Q 2 ) at Q 2 = 0, which defines the 
square of the neutron charge radius, is determined pre- 
cisely by the scattering of thermal neutrons from atomic 
electrons Q. 

This paper reports on new measurements of G^/G^ 
at low Q 2 in the vicinity of the maximum of G^, us- 
ing a longitudinally polarized electron beam incident on 
a vector-polarized 2 H target internal to the South Hall 
Ring at the MIT-Bates Linear Accelerator Center. The 
BLAST detector was used to detect quasielastically scat- 
tered electrons in coincidence with recoil neutrons over a 
range of Q 2 between 0.10 and 0.55 (GeV/c) 2 . 

The differential cross section for the 2 H(e, e'n) reaction 
with polarized beam and target can be written [13, EH El 
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where a unp is the unpolarized differential cross section, 
P z =n + — ?i_ and P zz = ?i + +n_ — 2?in are the vector and 
tensor polarizations of the deuteron target defined by the 
relative populations n m of the three deuteron magnetic 
substates with respect to the deuteron orientation axis, 
m = +1,0,-1, respectively, and P e is the longitudinal 
polarization of the electron beam. 

With BLAST, all of the polarization observables Ai 
in Eq. ^ have been measured for the first time with 
precision in a single experiment. The beam-target vec- 
tor polarization observable A^ d is particularly sensitive 
to the neutron form factor ratio G^/G 7 ^ [12j- In the 
Plane Wave Born Approximation (PWBA) and with the 
deuteron in a pure S-state, the asymmetry A^ d can be 
written analogously to elastic scattering from the free 
neutron as 
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where 9* and <jf are the target spin orientation angles 
with respect to the momentum transfer vector and a, b, 
and c are known kinematic factors. This asymmetry has 
the largest sensitivity to G E when the momentum trans- 
fer vector is perpendicular to the target polarization, i.e. 
9* = 90°. 

However, there are sizable corrections to the asym- 
metry in Eq. (|3|), mainly at low Q 2 where they are 
dominated by final state interactions (FSI). The relative 
contributions of meson exchange currents (MEC), iso- 
bar configurations (IC) and relativistic corrections (RC) 
become more significant as the momentum transfer in- 
creases (see Fig.[T]). Extracting G E must be done by com- 
parison with theoretical asymmetries that include these 
effects. 

The effects of FSI can be monitored with the other 
polarization observables in Eq. The asymmetries 

A e , A d , and Aj d all vanish in the Born approxima- 
tion due to parity and time reversal conservation and 
remain very small (below 1%) even in the presence of 
FSI. This permits these observables to be used to iden- 
tify any false asymmetries in the experiment. FSI gives 
a sizable contribution to the target tensor asymmetry 
A d , which is insensitive to G E and otherwise close to 
zero in the quasifree limit. Figure Q] displays a Monte 
Carlo simulation of the reaction mechanism effects on 
the asymmetries A^ d (upper panel) and A d (lower panel) 
as a function of Q 2 along with the measured values. 
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FIG. 1: Measured (solid blue points) and calculated beam- 
target vector polarization observable A^, d (upper panel) and 
tensor asymmetry (lower panel) for the 2 H(e, e'n)p reac- 
tion at 850 MeV, a target orientation of 6 d = 31.6° into the 
left sector of BLAST, and with neutrons detected in the right 
sector. The colored curves are Monte Carlo simulations based 
on the deuteron electrodisintegration model of Ref. [TTI ] (dot- 
ted magenta = PWBA, short-dashed green = PWBA+FSI, 
solid red = PWBA+FSI+MEC+IC+RC) using standard pa- 
rameterizations for the nucleon form factors (see text). In 
addition, the corresponding curves for G% = (dash-dotted 
red) and for elastic scattering from the free neutron (dashed 
black line) are shown. 



The calculations use the standard dipole form factor 
G D = (1 + Q70.71)" 2 for G P E , G^/m p , and G^/fi n , 
and 1.91t/(1 + 5.6t)G d for G% [H, where fi p = 2.79, 
[i n = —1.91, and r = Q 2 /(4m 2 l ). The good agreement of 
the measured tensor asymmetry with the full model 
supports the calculations of FSI for a reliable extraction 
of G E from the beam-target vector asymmetry A^ d at 
the percent level. 

On the other hand, the corrections at low Q 2 to A^ d 
measured in the 2 H(e, e'p)n reaction in quasifree kine- 
matics arc negligible [12[ , which allows for a precise deter- 
mination of the product of beam and target polarizations 
P e P z along with the proton form factor ratio G P E /G P M in 
this reaction channel 14 1 . 

The BLAST experiment was designed to carry out 
spin-dependent electron scattering from hydrogen [l5[ 
and light nuclei. Details on the experimental setup can 
be found in [l^]. The internal target consisted of an 
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atomic beam source (ABS) combined with an open-ended 
storage cell through which the stored electron beam 
passed continuously [l7j|- The ABS produced polarized 
monoatomic deuterium gas in the storage cell with nu- 
clear vector (V+: m=l; V— : m=— 1) and tensor (T— : 
m=0) polarization states. In addition, the helicity h of 
the electron beam was flipped every injection cycle. Lin- 
ear combinations of the six charge-normalized yields Yhm 
define all five polarization observables in Eq. ((2]). The 
experimental value of the beam-vector polarization ob- 



servable A^ d can be written as 
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where Y to t is the total yield obtained by summing up all 
six combinations hm. A modest magnetic holding field 
was applied to define the polarization angle 9d within the 
horizontal plane and to minimize the depolarization of 
target atoms. The variation of 9d was carefully mapped 
over the extent of the target cell. The average value 
of 6d was determined along with the tensor polarization 
Pzz by comparing the simultaneously measured tensor 
asymmetries in elastic scattering from tensor-polarized 
deuterium [lfl with those expected at low Q 2 based on 
a parameterization of previous data fl9] . 

The BLAST detector is a toroidal spectrometer (8 sec- 
tors) with the horizontal sectors instrumented with wire 
chambers, aerogel Cerenkov counters, thin plastic tim- 
ing scintillators, and thick plastic scintillator walls for 
neutron detection. With the target polarization vector 
pointing into the left sector, the neutron detection effi- 
ciency was augmented in the right sector covering the 
kinematic region most sensitive to the neutron form fac- 
tor ratio, as indicated by the sin#* term in Eq. ([3]). The 
detection of neutrons in the left sector was primarily used 
to independently verify the determination of P e P z from 
the 2 H(e, e'p)n reaction. 

The selection of (e, e'n) events is very clean; the num- 
ber of proton tracks misidentificd as neutrons is negligi- 
ble, due to the highly efficient charged particle veto pro- 
vided by the thin scintillator bars and the large-volume 
drift chambers in front of the neutron detectors. A set of 
cuts applied on the time correlation between the charged 
and the neutral track, and on kinematic constraints for 
the electrodisintegration process, was employed to iden- 
tify the quasielastic (e, e'n) events. The background from 
scattering off the aluminum target cell walls, measured 
with a hydrogen (empty) target, is less than 4% (3%) of 
the normalized yield obtained with deuterium. 

The corrected asymmetries were compared to Monte 
Carlo simulations based on the deuteron electrodisinte- 



gration model 111, for which events were generated ac- 
cording to the unpolarized cross section and weighted 
event- by-event with the spin-dependent terms in Eq. l[2"]). 
The acceptance-averaged asymmetry A^ d was simulated 



for different values of / G^ and compared to the ex- 
perimental values. In order to extract the best value of 
the form factor ratio for each Q 2 bin, a x 2 minimization 
was performed independently with respect to the missing 
momentum of the reaction and the angle of the neutron 
in the hadronic center-of-mass system. Both extractions 
produced consistent results. 

The data reported here were acquired in two separate 
runs in 2004 and 2005, corresponding to a target polar- 
ization angle of 31.64° ± 0.43° and 46.32° ± 0.45°, re- 
spectively. With a total accumulated beam charge of 451 
kC (503 kC) in the first (second) data set, final samples 
of 268,914 (205,252) coincident electron-neutron events 
were collected. The average product of beam and tar- 
get polarization determined from the 2 H(e, e'p) reaction 
was P e P z = 0.5796±0.0034(stat)±0.0034(sys) in the first 
and 0.5149±0.0043(stat)±0.0054(sys) in the second data 
set [14] . In comparison, the polarization product deter- 
mined from 2 H(e, e'n)p with neutrons detected in the left 
sector of BLAST corresponding to 9* « 0°, was found to 
be 0.587 ±0.019(stat) and 0.481 ± 0.026(stat) consistent 
with the above (e, e'p) results. The two data sets were 
treated as separate experiments producing two consistent 
results for the form factor ratio, which were combined for 
a final result. The data were divided into four Q 2 bins 
to determine G^/G^f with a comparable statistical sig- 
nificance (see Table HI . 

The systematic error of G^/G^ is dominated by the 
uncertainty of the target spin angle 9d- Other system- 
atic uncertainties include that of the beam-target polar- 
ization product P e P z , the accuracy of kinematic recon- 
struction, as well as the dependency on software cuts. 
The systematic uncertainties were evaluated individually 
for each Q 2 bin and data set by combining the errors 
from each source, taking covariances into account; the 
correlated and uncorrelated error categories of the two 
measurements were then combined for a resulting sys- 
tematic error of each bin. False asymmetries were stud- 
ied with the observables and A^ d and found to be 
consistent with zero. Radiative corrections to the asym- 
metries calculated in a PWBA formalism using the code 
MASCARAD [U are <1% and therefore also neglected. 
The uncertainties of the reaction mechanism and FSI cor- 
rections, which are small compared to the experimental 
errors, are not included in the systematic error. 



Q 2 /(GeV/c) 2 


(Q 2 )/(GeV/c) 2 




0.10-0.17 


0.142 


0.0505 ± 0.0072 ± 0.0031 


0.17-0.25 


0.203 


0.0695 ± 0.0084 ± 0.0039 


0.25-0.35 


0.291 


0.1022 ± 0.0127 ± 0.0046 


0.35-0.55 


0.415 


0.1171 ± 0.0182 ± 0.0052 



TABLE I: Results for the extracted neutron form factor ra- 
tio finG^ I 'G n M (fi n = Gm(0) = —1-91) with statistical and 
systematic errors, respectively. 
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The world's data on G^. from double-polarization ex- 
periments [f| are displayed in Fig.[2]along with the results 
of this work. All of the polarization data were experi- 
mentally determined as electric to magnetic form factor 



ratios. We used parameterization [2l| for G 7 ^, which 



is in good agreement with recent measurements |22(, to 



determine G^ from BLAST and to adjust the previously 
published values. The data from a variety of experiments 
are consistent and remove the large model uncertainty of 
previous Gg extractions from elastic electron-deuteron 
scattering [23J. The new distribution is also in agree- 
ment with G" extracted from the deuteron quadrupolc 
form factor [241 ] . 

The measured distribution of G% can be parameter- 
ized as a function of Q 2 based on the sum of two dipoles, 
EiOi/(l + Q 2 /h) 2 (i=l,2), shown as the BLAST fit in 
Fig. [2] (blue line) with a onc-sigma error band. With 
Ge(0) — and the slope at Q 2 = constrained by 
(r 2 ) = (-0.1148 ± 0.0035) fm 2 9], one parameter is 
fixed, resulting in a-y = —a 2 = 0.095 ± 0.018, bi = 
2.77 ±0.83, b 2 = 0.339 ±0.046 and cov(a 1 ,6 1 ) = -0.014, 
cov(ai,6 2 ) = 0.0008, cov(6i,6 2 ) = -0.036 with Q 2 in 
units of (GeV/c) 2 . The parameterization [25[ (magenta 
dash-dotted line) is based on the form introduced in [2lj 
with an additional bump structure around 0.2 — 0.4 
(GeV/c) 2 . Also shown are recent results based on vector 
meson dominance (VMD) and dispersion relations (red 
short-dashed 0] and green long-dashed lines Q ) , and the 
prediction of a light-front cloudy bag model with rela- 
tivistic constituent quarks [f| (cyan dotted line). 

The new data from BLAST do not show a bump struc- 
ture at low Q 2 as previously suggested 2l|, [25j. The 



BLAST data are in excellent agreement with VMD based 
models [3, HI and also agree with the meson-cloud calcu- 
lation [6(. The improved precision of the data at low 
Q 2 provides strong constraints on the theoretical under- 
standing of the nucleon's meson cloud. 

We thank the staff of the MIT-Bates Linear Accel- 
erator Center for delivering high quality electron beam 
and for their technical support, and A. Bernstein for sug- 
gesting the form of the BLAST fit. This work has been 
supported in part by the US Department of Energy and 
National Science Foundation. 
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